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Abstract 

Background: Calcium signals ([Ca^^j) direct many aspects of embryo development but their regulation is not well 
characterised. Ryanodine receptors (RyRs) are a family of intracellular Ca^^ release channels that control the flux of 
Ca^^ from internal stores into the cytosol. RyRs are primarily known for their role in excitation-contraction coupling 
in adult striated muscle and ryr gene mutations are implicated in several human diseases. Current evidence 
suggests that RyRs do not have a major role to play prior to organogenesis but regulate tissue differentiation. 

Findings: The sequences of the five zebrafish ryr genes were confirmed, their evolutionary relationship established 
and the primary sequences compared to other vertebrates, including humans. RyRs are differentially expressed in 
slow {ryrlo), fast (ryrS) and both types (ryrlb) of developing skeletal muscle. There are two ryr2 genes {ryr2o and 
ryr2b) which are expressed exclusively in developing CNS and cardiac tissue, respectively. In addition, ryrS and ryr2o 
mRNA is detectable in the initial stages of development, prior to embryonic axis formation. 

Conclusions: Our work reveals that zebrafish ryr genes are differentially expressed throughout the developing 
embryo from cleavage onwards. The data suggests that RyR-regulated Ca^^ signals are associated with several 
aspects of embryonic development, from organogenesis through to the differentiation of the musculoskeletal, 
cardiovascular and nervous system. These studies will facilitate further work to explore the developmental function 
of RyRs in each of these tissue types. 



Background 

Transient changes in the concentration of intracellular 
calcium ions ([Ca^'^]i) act as a powerful signal that is 
crucial for the establishment of form and function in 
the embryo. Detailed imaging studies have revealed that 
the spatial and temporal organisation of Ca^^ signals 
during embryogenesis are associated with many of the 
major phases of development, from early cell division to 
the differentiation of tissues. Despite their importance 
little is known about the generation and regulation of 
embryonic Ca^"^ signals. A comprehensive understanding 
of the pathways that regulate [Ca^"^]i during develop- 
ment is essential to understand the functional relevance 
of these signals in the embryo. 

Ryanodine receptors (RyR) are a family of intracellular 
Ca^"^ release channels that regulate the entry of Ca^"^ 
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into the cytosol from the intracellular organelles (the 
endoplasmic and sarcoplasmic reticulum). The RyR is a 
large homotetrameric protein (approximately 2,200 
kDa), each subunit is comprised of a large N-terminal 
cytoplasmic domain which modulates the gating of the 
channel, as well as luminal and transmembrane span- 
ning (TM) domains. In mammals, there are three dis- 
tinct ryr genes {ryrl, ryr2 and ryrS) that encode three 
differentially expressed RyR proteins. RyRl and RyR2 
are expressed predominantly in skeletal and cardiac 
muscle respectively, whilst RyR3 is found in many tis- 
sues at relatively low levels [1]. The primary role of the 
RyR is to increase [Ca^'^]i during excitation-contraction 
coupling (E-C coupling) in both skeletal and cardiac 
muscle. In humans, mutations in the ryrl and ryr2 
genes cause skeletal myopathies and cardiac disease 
respectively [2]. The generation of mouse knockout lines 
has provided some insight into the role of the receptors 
in the developing tissues of intact animals. Homozygous 
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mice from ryrr {skrrml) and ryr2~ knockout lines 
display gross morphological defects in either the skeletal 
muscle (ryrl) or heart tube (ryr2) and die at the perina- 
tal or embryonic day 10 (ElO) stages respectively [3,4]. 
In contrast, ryr^^"^"^ knockout mice appear to have no 
gross developmental defects and evidence suggests that 
RyR3 act to augment the [Ca^'^]i response of the other 
RyR isoforms in striated muscle [5,6]. The observation 
that RyR expression does not occur until relatively late 
in mammalian development [7], coupled with the fact 
that knock out lines are not lethal at very early stages 
has been interpreted as indicating that RyRs do not 
function during initial development. 

The zebrafish has been used extensively as a model 
for vertebrate development. The rapid development ex 
utero and embryonic transparency has proved advanta- 
geous for imaging the spatial and temporal organisation 
of Ca^"^ signals. These signals are involved in many of 
the early embryonic events; the initiation of fertilisation 
(0 hours post fertilisation, hpf), the early cell divisions 
associated with the cleavage period (up to 2 hpf) and 
the more extensive cellular rearrangements that occur 
in the blastula period (up to 5 hpf) (as reviewed in [8]). 
Evidence suggests that the release of Ca^^ from intracel- 
lular stores via the phosphatidylinositol (PI) signalling 
pathway is largely responsible for these early transient 
changes in [Ca^'^]i [9]. RyRs have not been implicated in 
very early developmental events and their expression 
during these stages has not been documented. Fluxes in 
the levels of embryonic [Ca^"^]i continue to occur 
throughout gastrulation (up to 10 hpf) (as reviewed in 
[8]). Initially changes in [Ca^"^]i occur as localised events 
but, as gastrulation progresses, co-ordinated waves of 
Ca^"^ signalling appear across the embryo. These later 
signals are proposed to coordinate a wide range of cel- 
lular movements (epiboly, involution, convergence and 
extension) that give rise to the embryonic body plan. 
Release of Ca^^ from intracellular stores via the PI sig- 
nalling pathway is again implicated at this stage; how- 
ever, the contribution to gastrulation from other Ca^^ 
signalling pathways remains undefined. Finally, the seg- 
mentation period (from 10 up to 24 hpf) is charac- 
terised by organogenesis and the emergence of the body 
systems. The Ca^"^ signals that occur during segmenta- 
tion are again more localised and typically associated 
with developing tissues. Transient changes in [Ca^"^]i 
have been recorded within the nervous system, somites 
and cardiac tissue [10-12]. Several studies in zebrafish 
have shown that inhibition of ryanodine receptor func- 
tion, using both pharmacological and genetic inhibitors, 
leads to impaired excitation-contraction coupling and 
gross morphological defects in the skeletal muscle, sug- 
gestive of a role in the development of this tissue 
[11,13,14]. 



This study set out to acquire a more comprehensive 
understanding of ryanodine receptor expression in early 
vertebrate development, using the zebrafish as an in vivo 
model. Our initial work confirmed the sequence of the 
five zebrafish ryr genes, established their evolutionary 
relationship to those in other vertebrate species and 
provided a direct comparison between the structural fea- 
tures of the primary protein sequences found in the zeb- 
rafish and mammals. An overview of ryr gene expression 
during zebrafish embryogenesis will inform work aimed 
at establishing the developmental significance of this 
family of Ca^^ -release channels. Therefore we conducted 
a comprehensive temporal and spatial analysis of ryr 
mRNA expression in the embryo using a combination 
of semi-quantitative PGR and wholemount in situ hybri- 
disation. We observed strong maternal expression of ryr 
mRNA [ryrS and ryr2d) during the cleavage and blastula 
periods suggestive of a novel role in early development. 
At 24 hours post fertilisation (hpf) ryr la, ryr lb and ryrS 
are expressed in skeletal muscle, whereas ryr2a is loca- 
lized to the central nervous system (CNS) and ryr2b is 
found exclusively in the cardiac muscle. Our study sug- 
gests that RyR channels have a role in early develop- 
ment prior to organogenesis as well as in the 
differentiation of different cell types. 

Methods 

Animal procedures 

Wildtype (WT) zebrafish strains {Tubingen and Tupfel 
long fin) were bred and raised in-house at the zebrafish 
facility of Queen Mary College, University of London, 
UK, as described previously [15]. Smoothened (smo) 
mutants were received as a gift from Prof. Simon 
Hughes (King's College London, UK). Embryos were 
collected by natural spawning and staged according to 
Kimmel and colleagues [16], given in the text as stan- 
dard developmental time at 28.5°C (hours post fertilisa- 
tion, hpf). Work on zebrafish embryos (prior to 
independent feeding) is exempt under the U.K. Animals 
(Scientific Procedures) Act 1986 and does not require 
ethical approval. 

Genomic analysis and gene prediction of ryr genes 

To identify zebrafish ryr genes encoding RyRs as 
described in Hirata and colleagues [13], protein 
sequences corresponding to the human RyR family 
[http://www.ensembl.org/Homosapiens/familyview7famil- 
y=ENSF00000000736] were used as template for search- 
ing in GenBank and the zebrafish genome resources at 
Ensembl http://www.ensembl.org/Danio_rerio/; Zv7 
Ensembl assembly using a protein tBLASTn approach. 
Several hits were identified, which were then compared 
to the 3'-UTR sequence of the EST to identif)^ the corre- 
sponding genomic region. The gene predictions 



Wu et al. BMC Research Notes 201 1, 4:541 
http://www.biomedcentral.eom/1756-0500/4/541 



Page 3 of 1 7 



program, Genwise (http://www.ebi.ac.uk/Wise2/) was 
applied to reveal genomic contigs and full sequence 
relating to the ryrla, ryr2a, ryr2b and ryrS in zebrafish 
(Table 1). For ryrla and ryrS, based on the tBLASTn 
result, the Accession numbers (Table 2) and annotated 
gene sequences for these isoforms were selected. For 
ryr2ay owing to the limitations to the resources available 
on the database, 5 contigs (NA_1034, NA_1216, 
NA_3083, NA_1397 and NA_1713) were identified with 
nucleotide similarity to the human ryr2, based on the 
approach described above. Bioinformatics database 
search approaches using human ryr2 did not reveal any 
fully annotated ryr2b gene sequences in zebrafish. How- 
ever, database blast results located a gene (EntrezGene 
Name LOC568506 on chromosome 17:18.7 m) encoding 
805 residues. The partial ryr2b gene located in a 19.98 
kb region in the Ensembl Zv7 assembly shared very lim- 
ited conserved symmetry with the human ryr2. An 
alignment approach revealed that the zebrafish ryr2b 
gene was reverse transcribed and the region of 19.98 kb 
encoding 805 residues at the C-terminus end shows 
high peptide identity to the human RyR2 C-terminus. 
To help to find the missing large N-terminal cytoplas- 
mic domain of the RyR2b peptide, Genscan features on 
Ensembl Zv7 assembly were employed to reveal an extra 
1014 residues, giving a total of 1819 residues from the 
C-terminal end of the peptide. In order to search further 
for the missing portion of the RyR2 peptide towards the 
N-terminus, a genomic region of 107,783 bp spanning 
the C-terminal end of ryr2b gene to the start of neigh- 
bouring actn2 gene [Ensembl: ENSDARG00000071090], 
was reverse complemented using REQSEV program 
http://bioweb.pasteur.fr/seqanal/interfaces/revseq.html 
before performing a pairwise BLAST alignment with 
Gasterosteus aculeatus (stickleback) RyR2 peptide 
sequence on NCBI (http://www.ncbi.nlm.nih.gov/blast/ 
bl2seq/wblast2.cgi) to confirm its sequence identity. 
Based on the alignment result, the full length zebrafish 
RyR2b peptide was first constructed with reference to 
Takifugu rubripes (pufferfish) and stickleback RyR2 pep- 
tide sequences using the Genewise program (http:// 



www.ebi.ac.uk/Wise2/). The complete RyR2b peptide 
was finalised by the alignment of the translated genomic 
region using SHOWORF program http://bioweb.pasteur. 
fr/seqanal/interfaces/showorf.html. The sequence for the 
ryr2b gene was deposited in EMBL Nucleotide Sequence 
Database (Accession no. FR822741). For ryrlby sequence 
was obtained via the NCBI CoreNucleotide search web- 
site [NCBI - AB247454] and isoform-specific PCR pri- 
mers (set 3) were used, as described previously [13]. 
Isoform-specific PCR primers for ryrla, ryr2a, ryr2b and 
ryr3 were designed by retrieving the nucleotide 
sequences based on divergent regions of the RyRs iden- 
tified from peptide sequence (Figure 1). These regions 
were re-confirmed by both BLAST against different 
databases (NCBI nucleotide blast, Ensembl and ZFIN) 
and DNA sequencing. 

Phylogenetic Analysis 

Protein sequences equivalent to the RIH_assoc 
(pfam08454) and RR_TM4-6 (pfam06459) domains in 
the five zebrafish RyR protein sequences (RyRla, RyRlb, 
partial RyR2a, RyR2b and RyR3) and other well charac- 
terised vertebrates RyR homologue sequences were 
extracted from Ensembl and GenBank databases (Table 
2). The RyR sequences were pre-aligned using the Clus- 
talW alignment program to remove any gaps generated 
within the sequences, followed by a multiple RyR 
sequence alignment using the T-coffee program avail- 
able from the European Bioinformatics Institute web- 
page (http://www.ebi.ac.uk/Tools/t-coffee/index.html) 
with the default parameters. The output from the multi- 
ple sequence alignment result obtained from T-coffee 
was used subsequently as the template for the genera- 
tion of a Guide Tree using the equal angle method 
available on the SplitsTree4 program with its default 
parameters (http://www.splitstree.org/for free download). 

Semi-quantitative end-point PCR 

Total RNA was isolated from nine stages of zebrafish 
development: 4- to 128-cell stage (1-2.25 hpf), 50%-epi- 
boly stage (5.3 hpf), 100%-epiboly stage (10 hpf), 3- to 



Table 1 Summary of the ryr genes identified in zebrafish and their identity to the human orthologues 



Genes 


cDNA length (bp) 


Numbers of Exons 


Genome Location 


Coding Region (aa) 


Identity to Human isoforms# 


ryrla 


15,195 


108 


Chromosome 10: 29.76 m 


5,064 


76% 


ryr lb 


15,231 


100* 


Chromosome 18: 36.78 m 


5,076 


76% 


ryr2a 


N/A 


N/A 


5 contigs so far; scaffolds 


2,829 


86% 


ryr2b 


14,748 


121* 


Chromosome 17: 18.70 m 


4,916 


63% 


ryrS 


15,122 


98 


Chromosome 20: 38.1 1 m 


4,863 


77% 



The zebrafish /yr/a, ryrib, ryr2b and ryrS genes are composed of 108, 100, 121 and 98 exons, and encode proteins of 5064, 5076, 4916 and 4870 amino acids 
(aa), respectively (Data from Zv7 Ensembl and BLAT Search). To date, five genomic contigs were identified for ryr2a, which spans 2829 amino acids, based on a 
protein tBLASTn searching approach using human RyR2 peptide 
* = BLAT RESULTS and # = ClustalW Alignment Score 
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Table 2 Sources of ryanodine receptor sequences used for phylogenetic analysis 



RyRs 


Species 


Transcripts 


Peptides 


RyRI 


Chimpanzee 


LINbr 1 K 1 00000020262 


ITN 1 c n~rn nr\r\r\r\r\r\ 1 o~7/i i 

LINbP 1 KPOOOOOOl 8/41 


RyRl 


Human 


LNb 100000359596 


LI\lbP00000352608 


RyRI 


Mouse 


LNbMUb 1 0000003281 3 


LN3MU3P0000003281 3 


RyRl 


Pig 


N 1 n r\r\ i r\r\ i n /i 

NP_001001534 




RyRI 


Rabbit 


PI 1 716 




KyKl_(l ) 


Medaka 


rKiC/^ni T r\r\r\r\r\r\r\ci r\r\^ 

LNbUKL 1 00000008002 


LN3UKLP00000008001 


KyKl_(2) 




FNiC/^ni Tr\r\r\r\r\r\r\^ '->r\r 

LNbUKL 1 00000001 305 


LNbUKLPOOOOOOOl 304 


KyKl_(l j 


ruttertisn 


LI\Jb 1 KU 1 000000391 20 


iTNic — rni \r>r\r\r\r\r\r\'D or\or\ 

LINb 1 KUP00000038980 


KyKl_(z) 




LNb 1 KU 1 000000435/1 


LINb 1 KUP00000043426 


DwD 1 M ^ 

KyK1_(l j 


Sticl<lebacl< 


LN bG AL 1 000000 1/168 


LINbGALPOOOOOOl /1 34 


KyKl_(2) 




LN3GAC 1 0000002/379 


LN5GAGP00000027327 


RyRI a 


Zebrafish 


LN3UAK 1 0000003601 5 


LN3UAKP00000032856 


RyRl b 




Ad247454 




RyR2 


Cliicl<en 


LINbbAL 1 0000001 7582 


LINbGALPOOOOOOl 756 1 


RyR2 


Cliimpanzee 


LINbP 1 K 1 00000003929 


LINbP 1 KP00000003630 


RyR2 


Dog 


ENSCAFT00000015962 


ENSCAFP00000014764 


RyR2 


Human 


ENST00000S60064 


ENSP00000S5S174 


RyR2 


Macaque 


ENSMMUT00000001512 


ENSMMUP0000000142S 


RyR2 


Mouse 


ENSMUST00000021750 


ENSMUSP00000021750 


RyR2 


Rabbit 


PS0957 




RyR2 


ruttertisli 


GeneScan 




RyR2 


Sticl<lebacl< 


ENSGACT00000020252 


ENSGACP0000002021S 


RyR2b 


Zebrafisli 


EMBL Nucleotide Sequence Database (Accession no. FR822741) 


RyR3 


Cliicl<en 


ENSGALT00000005S52 


ENSGALP00000005S42 


RyR3 


Cliimpanzee 


ENSPTRT00000049180 


ENSPTRP0000004192S 


RyR3 


Dog 


ENSCAFT00000012974 


ENSCAFP00000012007 


RyRS 


Horse 


ENSECAT00000016127 


ENSECAP0000001S008 


RyR3 


Human 


ENST00000S62047 


ENSP00000S54852 


R\/R^ 


ICO 
IVIUUbc 


ENSMUST00000091818 


ENSMUSP00000089426 


RyRSJl) 


Medal<a 


ENSORLT00000022S71 


ENSORLP00000022S70 


RyR3_(2) 




ENSORLT000000211SS 


ENSORLP000000211S2 


RyRSJl) 


Pufferfish 


ENSTRUT00000046S40 


ENSTRUP00000046185 


RyR3_(2) 




ENSTRUT000000S6458 


ENSTRUP000000S6S27 


RyRSJl) 


Sticl<lebacl< 


ENSGACT00000008056 


ENSGACP000000080S7 


RyR3_(2) 




ENSGACT00000010590 


ENSGACP00000010568 


RyRS 


Zebrafisli 


ENSDART0000004655S 


ENSDARP00000046552 



6-somite stage (11-12 hpf), 18-somite stage (18 hpf), 24 
hpf, 72 hpf, 144 hpf stages and a 6-month old adult, 
using TRIZOL® Reagent (Invitrogen). Total RNA (5 
was subjected to 30 minutes treatment with DNase I 
(Promega) at 37°C prior to cDNA synthesis performed 
using a Cloned AMV First-Strand cDNA Synthesis kit 
(Invitrogen) with random hexamers. For semi-quantita- 
tive end-point PGR, 1.5 (ig of cDNA templates made 
from the developmental series were used alongside with 
the five sets of isoform-specific oligonucleotide primer 



pairs for amplification and normalised against the 
expression of fi-actin (Table 3). 

Whole mount in situ hybridisation 

In situ hybridisation of whole-mounted zebrafish 
embryos was performed as described previously [17]. 
Briefly, ryrla, ryrlb, ryr2a, ryr2b and ryrS sense and 
anti-sense digoxygenin (DIG) labelled RNA probes cov- 
ering 1,405 bp, 1,361 bp, 1,313 bp, 1,266 bp and 1,549 
bp, were generated (Table 3). Briefly recombinant vector 
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Human RyRl 












-S50JS 


Human RyR2 














- N4967 


Human RyR3 






f — 1 








-C487(( 


Zebrsifish Ry R L a 














-05064 


Zebrafish RyRlb 














-GS076 


Zebrafish R\ Rhi 












M ■ ■ Mt 


— Nxxxx 


Zebrafisb RyRZb 














-S4!>16 


Zebrafish RyR3 














-G4)i(i3 


KEY: 



■I MIR Doniiiin ■ RyR Donutin fl EF-tland (Ca^^ BitidiiiK> 

m RIH Domain ■ RIH A^KOi^iati-d UZ RvR TM -*-r» & (Pore formmj:) 

■ SPRY Dumjtin M Sei|ucnce Not Avails bte 500 ;immri ^ictilji 



(B) 



MS 

OLVMTVGLLAVVVVLYTVVAFN 
OLVMTVGLLAVVVYLYTVVAFN 
OLVMIVGLLAVVVVLYTVVAFN 

RyHI Pig OLVMIVGLLAVWYLYTWAFN 

Ryfil (1) Kedaka. ...... OLMMT\/GLLAVWYLYTVVAFN 

RyR\ (2) Medska y,my.y!.>!.MU)i.y,K)imx^xu 



Ryfll Chimpanzees, 

RyRl Human 

RyRl Mouse 



Ryfll(l) Pufferftsh. 

mH2) Pufferftsh. 

Ry«1(1) Stickleback. 

Ry«1C2) Stickleback. 



OLUMIVGLLAVVVYLYTVVAFN 
OLW«rVGLLAVVVYLYTVVAFN 
OLMMTVGLLAVVVYLYTVVAFN 
OLiflWTLGLLAVVVYLYTVVAFI 



RyRl a Zebrafish.,, OLWHrVGLLAVVVYLYTVVAFW 

RyRIb Zebrafish. ...... OLLMTVGLLAVVVYLYTVVAFN 

RyR2 Chicken OLVLTVGLLAVmLYTVVAFN 

RyR2 Chimpanzees OLVLTVGLLAVVVYLYTVVAFN 

RyR2 Dog. {JLVLTVGLLAVVVYLYTVVAFN 

RyH2 Hyman OLVLTVGLLAVVVYLYTVVAFN 

RyR2 Macaque OLVLTVGLLAVVVYLYTVVAFH 

RyR2 Mouse OLVLTVGLLAVVVYLYTVVAFN 

RyR2 Rabbit OLVLTVGLLAVVVYLYTVVAFN 

RyR2 Puff erf i sh. ..... . OLVLTVGLLAVVVYLYTVVAFN 

RyR2 St I Ck I eback. ..... OLVLTVGLLAVVVYLYTVVAFN 

RyR2a Zebr^. i . OLVLTVGLLAVVVYLYTVVAFN 

RyR2b Zebrafish OLVLTVGLLAVVVYLYTVVAFN 



Pore Mm 

GGGIGD RVVFDITFFFVIVILLAIiQGLlI 

GGG I GD RVVFD I TFFFV E V I LLA 1 1 QGLl I 

GGGIGD RVVFDITFFFFViVILLAIIOGLII 

GGGIGD RVVFDiTFFFVIVILLAIIQGLII 

GGGIGD RVVFDITFFFVIVJLLAMQGLII 

GGGIGD RVVFDITFFFFVIVILLAIIOGLII 

GGGIGD RVVFDITFFFFVIVILLAIIOGLII 

GGGI GD RVVFD I TFFFFV I V I LLA 1 1 OGLl I 

GGGIGD RVVFDITFFFViVILLAIIOGLlI 

GGGIGD RVVFDITFFFFVIVILLAIIOGLII 

GGGIGD RVVFDITFFFFVEVILLAIIQGLTI 

GGGIGD RVVFDITFFFFVIVILLAIIQGLII 

GGGIGD RIIFDITFFFFVIVILLAIIQGLII 

GGGIGD RIIFDITFFFVIVILLAIIQGLII 

GGGIGD RllFDITFFFFVfVILLAIIQGLH 

GGGIGD REIFDITFFFVIVILLAIIQGLIl 

GGGIGD R! IFDITFFFVIVILLAI I QGLl I 

GGGIGD Rl IFDITFFFVIVILLAI lOGLH 

GGGIGD RIIFDITFFFVIVILLAIIQGLII 

GGGIGD RMIFDITFFFVIVILLAIIQGLII 

GGGIGD RIIFDITFFFVIVILLAIIQGLII 

GGGIGD RIIFDITFFFVfVlLLAnOGLlI 

GGGIGD RIIFDITFFFFVIVILLAIIQGLII 



RyR3 Bat 

RvR3 Ghioken. 

RyR3 Ch impanzees. . . . 

Ryfl3 Dog. 

RyR3 Horse 

RyR3 Human. . , 

RyR3 Mouse, . . , 

RyR3 Rabbit 

RyH3C1) Hedaka 

Ryll3(2) Medaka 

Rya3Cl) Pufferfish.. 
Ryl?3(2) Pufferfish.. 
Ryl?3Cl) Stickleback. 
RyR3(2) Stickleback. 
RyR3 Zebrafish. 



OLVLTVGLLAVVVYLYTVVAFN GGGIGD RIVFDITFFFVIVILLAI lOGLJ I 

OLVLTVGLLAVVVYLYTVVAFN GGGIGD RIVFDITFFFVIVILLAI IQGLlI 

OLVLTVGLLAVVVYLYTVVAFN GGG I GD RI VFD I TFFFFV I V I LLA 1 1 QGLl I 

OLVLTVGLLAVVVYLYTVVAFN GGGIGD RIVFDITFFFVIVILLAI IQGLl I 

OLVLTVGLLAVVVYLYTVVAFN GGGIGD RIVFDITFFFVIVILLAI lOGLl I 

OLVLTVGLLAVVVYLYTVVAFN GGGIGD RIVFDITFFFVtVILLAHQGLl I 

OLVLTVGLLAVVVYLYTVVAFN GGGIGD RIVFDITFFFVIVILLAIIQGLl I 

OLVLTVGLLAVVVYLYTVVAFN GGG I GD R I VFD I TFFFV I V I LLA 1 1 QGL 1 1 

OLVLTVGLLAVVVYLYSVVAFN GGGIGD RILFDITFFFVIVILLAI IQGLl I 

OLVLTVGLLAVVVYLYTVVAFN GGG I GD R I LFD I TFFFFV E V I LLA 1 1 QGL 1 1 

OLVLTVGLLAVVVYLYTVVAFN GGGIGD RILFDITFFFVIVILLAI IQGLH 

OLVLTVGLLAVVVYLYTVVAFN GGGIGD RILFDITFFFVIVILLAI I QGLl 1 

OLVLTVGLLAVVVYLYTVVAFN GGGIGD RILFDITFFFVIVILLAI IQGLl I 

OLM.TVGLLAVVVYLYTVVAFN GGGIGD RILFDITFFFVIVILLAI lOGLl I 

OLVLTVGLLAVVVYLYTVVAFN GGGIGD RILFDITFFFVIVILLAI JQGLI I 

V V 

A B 



mfftiTCTit 



1 



0 



mil 



MR 



iiiiii: 



Tmin! 



Port? 



mnmiiiiiiiiiiiiii 



MIA 



Figure 1 Analysis of the primary sequence of the zebrafish RyRs. (a) Sciiematic representation of tine protein sequence of zebrafisli RyRs 
compared to the iiuman isoforms. The figure siiows the tiiree iiuman RyRs at tine top and tine five zebrafisli RyRs below. Sequence equivalent 
to MIR and RIH domains identified in all five zebrafisli RyRs. Three SPRY domains, four copies of tine RyR repeated domain, a single copy of tiie 
RIH associated domain and a putative EF liand domain were also identified witliin the zebrafisli RyRs. Transmembrane domains (TM 4-6) and 
putative pore forming regions (TM 8-10) were also present, (b) Alignment of 36 RyR ortiiologues to siiow the sequences of the pore and 
transmembrane regions (M8, pore, MIO). Letters A and B indicate regions not shown in alignment. Bar graph indicates the number of different 
amino acids among the different isoforms at the indicated residues. 



Wu et al. BMC Research Notes 201 1, 4:541 
http://www.biomedcentral.eom/1756-0500/4/541 



Page 6 of 1 7 



Table 3 Primer sequences used for PCR and WISH protocols 



Genes 


(A) PCR Primers (5 to 3 ) (B) WISH Primers (5 to 3 ) 


ryr lo F 


\ LA 1 GLLL 1 LAbLA 1 GL 


fyflo R 


1 1 GbAA 1 L 1 LLGLL 1 LLbAG 1 A 


lyi 1 u r 


rrTrrTrrTi^riTT AT(^r a cc ATr;TT ATr, a r, 


ryrlb R 


TCAATGTGGTOTGATCAGTGGGAACGGGA 


ryr2a F 


GAAGATGGCGAGAAGAAACCTG ATGTCTGACGCCGGTGAACGAT 


ryr2o R 


TGCAGCATCATOTOTCCGTO CAGCACAAAGTGAACAGTGACGG 


ryr2b F 


TCTAAATTCAGCCACCATGATGTGAGC 


ryr2b R 


TGTOCTCTCATGTACTCGCCC 


ryrS F 


AAGGAAAAGAATGAGTCGGGC 


ryrS R 


CACGCAICI 1 1 ICI 1 IGGGA 



templates (pGEM-T Easy-ryrs) were linearised with an 
appropriate restriction enzyme and then subjected to 
phenol/chloroform purification and alcohol precipita- 
tion. Purified, linear DNA (0.5 (ig) was used for the in 
vitro transcription reaction (a 20 (il reaction contained 
IX Transcription Buffer (Invitrogen), 2 (il DTT (0.1 M), 
2 (il nucleotide mix (1 mM GTP, 1 mM ATP, 1 mM 
CTP, 0.65 mM UTP and 0.35 mM DIG-ll-UTP), 50U 
placental ribonuclease inhibitor and lOU of RNA poly- 
merase). The mixture was incubated at 37°C for 4 hours 
before the removal of the original DNA template by 
incubating with 2U of DNase at 37°C for 1 hour. DIG- 
labelled cardiac myosin light chain 2 (cm/c2; [18]), myo- 
genic differentiation {myoD;[19]), nkx2.5 ([20]) and 
fluorescein-ll-UTP labelled myosin heavy chain 1 
{myhzl;[21]) anti-sense probes were also prepared under 
the same approach and used as positive markers. 

Zebrafish embryos were fixed in 4% paraformaldehyde 
(PFA)/phosphate buffered Triton XlOO (PBT) at 4°C 
overnight or at room temperature (RT) for 2 hours, 
then washed in 100% methanol and incubated at -20°C 
for at least 30 minutes. Embryos were rehydrated in a 
methanol gradient and treated with 10 (ig/ml of Protei- 
nase K/PBT for 1 to 20 minutes at RT depending on 
their stages prior to fixation in 4% PFA/PBT for 20 min- 
utes at RT (embryos before or at 50% epiboly stage were 
not treated with Protease K or re-fixed). Embryos were 
washed with PBT and pre-hybridised in hybridisation 
mix (50% formamide, 5X SSC (0.75 M NaCl and 75 
mM triNa citrate) at pH 5.0, 500 (ig/ml yeast RNA, 50 
(ig/ml heparin, 0.1% Tween20) for 2 hours at 65°C. The 
mixture was replaced by fresh hybridisation mix con- 
taining DIG-labelled or Fluorescein-labelled RNA probe 
at the appropriate dilution and incubated at 65°C over- 
night. The post hybridisation washes were carried out at 
65°C. The excess probe was removed by sequential 
washing in 25% formamide in 2X SSC, 2X SSC, 0.2X 
SSC at 65°C and finally in PBT at RT. Embryos were 
incubated in maleic acid buffer (0.1 M Maleic acid pH 



7.5 and 0.15 M NaCl) containing 2% Blocking Reagent 
(Roche; MAB) at RT for a minimum of 1 hour. Staining 
was carried out by replacing the MAB solution with the 
1:5000 MAB diluted alkaline phosphatise (AP) conju- 
gated anti-DIG antibody (Roche) or 1:4000 MAB diluted 
AP conjugated anti-Fluorescein antibody (Roche) at 4°C 
overnight. The embryos were washed with PBT and 
then detection buffer (0.1 M Tris-HCl pH 9.5, 0.1 M 
NaCl and 50 mM MgCl2 and 0.1% Tween20) at RT. 
The AP substrates used for colour development were 
BM purple (Roche), Fast Red tablets (Roche) and NBT/ 
BCIP (Roche). The colour reaction was terminated by 
washing in PBT and embryos were re-fixed in 4% PEA/ 
PBT for 2 hours at RT. 

Immunocytochemistry 

Wholemount immunocytochemistry (WICC) was con- 
ducted as described previously [22]. Briefly, embryos 
were fixed overnight at 4°C and all subsequent steps 
were performed at RT. Embryos were incubated in 
blocking buffer prior to incubation in primary antibody 
in 1% goat serum in phosphate buffered saline supple- 
mented with 0.8% triton XlOO (PBST) at room tempera- 
ture. The primary antibodies used were 34 C, F59 and 
MF20 at dilutions of 1:250, 1:10 and 1:100, respectively. 
Embryos were rinsed in phosphate buffer with 0.8% tri- 
ton (PBT) and incubated in a 1:1000 dilution of goat 
anti-mouse IgG Cy^^-5 linked secondary antibody 
(Amersham) made up in 1% goat serum in PBT over- 
night. Embryos were rinsed in PBT and stored in 50% 
glycerol/50% PBS. 

For the double immuno-labelling (WISH/WICC), the 
WISH labelled embryo was initially developed using 
reagents as described above, followed by antibody label- 
ling using a 1:10 dilution of F59 and a 1:5000 dilution of 
goat anti-mouse horseradish peroxidase (HRP) conju- 
gated secondary antibody (Merck) carried out as 
described in the Zebrafish Book [23]. Colour develop- 
ment by HRP was activated by the addition of 0.01% v/v 
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DAB and 6% v/v H2O2. Reaction was terminated by the 
addition of PBS and fixed overnight in 4% PFA at 4°C. 
For double-fluorescent labelling, the WISH labelled 
embryo was initially developed using Fast Red, followed 
by antibody labelling as described above. 

Embryos were mounted in 100% glycerol and for flat 
mounting the yolk of the stained embryo was initially 
removed. Brightfield and fluorescent single or Z-stack 
images were collected using a XIO, X20, X63 oil immer- 
sion or XlOO water immersion objectives on a Zeiss 
LSM 510 microscope and LSM Image Examiner soft- 
ware. LSM Image Browser (Version 4,2,0,121) software 
was used for image post-processing such as preparing 
three dimensional (3D) projected cross-sections from 
the acquired Z-stacks of images. For the observation of 
the Cy™5, Fast Red and Alexa Fluor® (Alexa Fluor® 
488 nm) staining, argon lasers 633 nm, 543 nm and 488 
nm lines and LP 650, LP 560 and BP 505-550 emission 
filter(s) were used, respectively. 

Results 

The relationship between ryanodine receptor genes in 
zebrafish and those of other vertebrates 

We conducted a tBLASTn search of the zebrafish geno- 
mic database (Zv7 Ensembl), using the protein 



sequences for human RyRl, RyR2 and RyR3 as queries, 
and detected four zebrafish RyRs homologues: ryrla, 
ryrlb, ryr2b and ryrS (as summarised in Table 1). ryr2a 
was identified from scaffold sequences by blasting the 
human RyR2 sequence in the updated NCBI database. 
Despite the fact that the sequence obtained from ryr2a 
is not a completed version based on the available infor- 
mation in current databases, the five genomic contigs 
(NA_1034, NA_1216, NA_3083, NA_1397 and 
NA_1713) identified to date show high similarity to 
other vertebrate RyR orthologues. Zebrafish ryr gene 
sequences display significant homology to their respec- 
tive human isoforms (Table 1). 

The evolutionary relationship of vertebrate RyR genes, 
including the five zebrafish sequences was examined 
(Figure 2). Three major monophyletic clades were iden- 
tified that correspond to RyRl, RyR2 and RyR3 subfami- 
lies. We found that the zebrafish RyRla, RyRlb, RyR2a, 
RyR2b and RyR3 all clustered with their respective gene 
families into the three monophyletic clades [24]. The 
zebrafish RyRla and RyRlb cluster to their respective 
subgroups with those of the other teleost species, Gas- 
terosteus aculeatus (stickleback), Oryzias latipes 
(medaka) and Takifugu rubripes (pufferfish). We identi- 
fied two ryr2 genes in zebrafish, which currently appears 




Human IPjR(J) 
Human IPjR(2) 
Human IPjR(l) 



Figure 2 Phylogenetic tree of the vertebrate RyR gene family. The evolutionary relationship between the zebrafish ryr orthologues and 
cloned full-length ryr genes of other vertebrate species was determined. RyR protein sequences were aligned using the T-coffee multiple 
sequence alignment. The phylogenetic tree was generated using the equal angled method available on the SplitTree4 software. The branch 
length of the tree is proportional to evolutionary distance. The family of human inositol 1,4,5-trisphosphate receptors (IP3R) was used as an out- 
group. Three monophyletic clades corresponding to (1) RyRl, (2) RyR2 and (3) RyR3 gene subfamilies were identified. A full description of the 
sequences used can be found in Table 2. through to adulthood. 
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to be the only species in which this gene duplication 
event can be observed. The zebrafish ryr2a gene appears 
to be more closely related to other teleosts (i.e. stickle- 
back and pufferfish) whereas ryr2b gene is more 
diverged. Although most teleosts appear to have two 
ryrS genes, we confirmed previous studies which report 
a single zebrafish ryrS gene [11,13,24]. 

Zebrafish RyRs contain many of the conserved structural 
domains with similarities to other vertebrates 

The molecular structure of the RyR protein has been 
explored extensively, partly in the drive to understand 
its regulation. Each monomeric RyR protein is approxi- 
mately 5,000 amino acids in length with a molecular 
mass of 565 kDa. The receptor has a large N-terminal 
cytoplasmic domain containing many regulatory binding 
sites (as reviewed in [25,26]) that modulate the gating of 
the channel pore located in the C-terminus. The N-ter- 
minus cytoplamic domain of RyR interacts with a host 
of regulatory proteins, such as calstabin and calmodulin. 
Physiological modulators of RyR function include ATP, 
Ca^^, Mg^^, cyclic ADP ribose, posttranslational modifi- 
cations (e.g. phosphorylation, oxidation) and pharmaco- 
logical substances (e.g. ryanodine, caffeine) [27]. 

We conducted an analysis to explore whether the zeb- 
rafish RyR protein sequences contain characterised con- 
served domains that may contribute to receptor 
regulation and their comparison to their human coun- 
terparts was analysed (Figure la). MIR (Mannosyltrans- 
ferase. Inositol 1,4,5-trisphophate receptor (IP3R) and 
RyR [pfam02815]) and RIH (RyR and IP3R Homology 
[pfam01365]) domains were identified within the N- 
terminal of all five zebrafish RyRs. The MIR and RIH 
domains are common to all the members of the intra- 
cellular Ca^"^-release channel super family [28,29]. The 
MIR domain has been suggested to have a ligand trans- 
ferase function and the RIH domain may form a binding 
site for IP3; however, very little is known regarding their 
role in receptor regulation to date [29]. All of the zebra- 
fish RyRs contain three SPRY domains (SPla and the 
RyR [pfam00622]), which have been proposed to interact 
with voltage gated channels [30,31]. Each zebrafish RyR 
also contains four copies of the RyR domains (RyR 
repeated domain [pfam02026]), a sequence unique to 
these channels [29]. Furthermore, zebrafish RyRs con- 
tain the eukaryotic RIH associated (RyR and IP3R 
Homology associated [pfam08454]) domain, which cur- 
rently has no known function. EF-hand motifs may have 
a functional significance in activation of the channel by 
Ca^"^ themselves. Putative Ca^"^ binding sequences, EFl 
and EF2 [32], have been identified in the receptor and 
are thought to be the major Ca^^ regulatory sites. We 
identified similar putative EF-hand motifs towards the 
C-terminal of the receptor and these appear to be 



conserved in all zebrafish and human isoforms. The RyR 
pore determines the conductance and ion selectivity of 
the channel; however, the structure of the TM and pore 
forming region is still unresolved. There are eight pro- 
posed TM sequences, of which the last six are suggested 
to form the Ca^"^ release channel [33]. According to the 
topological model by Du and colleagues six to eight TM 
sequences (i.e. M4a/M4b, M5, M6, M7a/M7b, M8 and 
MIO) were identified, with the M9 (pore segment) 
inserted between the M8 and MIO TM segments. This 
M9 region is proposed to act as the selectivity filter 
allowing Ca^"^ to transverse the membrane [34]. All five 
zebrafish RyRs contain RR TM 4-6 domains (RyR TM 
4-6 [pfam06459]). Comparison of sequence in the core 
pore-forming region demonstrate that all five zebrafish 
RyR sequences are highly conserved with other verte- 
brate species, although it should be noted that there are 
subtle changes at the single amino acid level (Figure lb). 

Zebrafish ryr mRNA expression can be detected from the 
earliest stages of development through to adulthood 

We examined the temporal expression of the ryr genes 
during nine key stages of zebrafish development using 
semi-quantitative end-point PGR (Figure 3a). Expression 
of ryr la was detected from 11-12 hpf (3- to 6-somite 
stage) through to adulthood, whereas ryr lb expression 
was evident from 18 hpf (18-somite stage) onwards. 
Results also showed that levels of mRNA expression for 
both ryr la and ryr lb were at the highest in 144 hpf and 
adult stages. The ryr2a gene showed strong maternal 
expression at 1-2 hpf (4- to 128-cell stage) and zygotic 
expression was observed from 18 hpf through to adult- 
hood. Expression of ryr2b mRNA was detected from 18 
hpf through to adulthood with relatively weak expres- 
sion observed at 24 hpf and 144 hpf stages. Finally, ryrS 
maternal expression was observed at 1-2 hpf, followed 
by weak zygotic mRNA expression from 5.3 hpf (50% 
epiboly stage) through to 18 hpf (18-somite stage), after 
which the levels became much stronger ryrla, ryrlb 
and ryrS are expressed exclusively in skeletal muscle 
during segmentation period 

We examined the spatial and temporal distribution 
pattern of the zebrafish ryr mRNA in the developing 
embryo from 12 to 24 hpf of development using WISH. 
The expression of ryr la, ryrlb and ryrS in the develop- 
ing somites was compared to the myogenic determinant 
(myoD), one of the earliest markers of myogenic com- 
mitment with a key role in regulating muscle differentia- 
tion (Figure 4a) [19]. 

We observed low levels of ryrla expression in the 
adaxial cells located on either side of the notochord 
beginning at 11 hpf (Figure 4a). Our findings demon- 
strate that in slow muscle ryrla appears at 11 hpf, 1 to 
2 hpf prior to fast muscle ryrlb mRNA expression. 
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(A) 



Genes 
ryrla 
ryrlh 
ryrla 
ryr2b 
lyrS 
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(B) 



Tissues 

BR BO Sizes (bp) 

279 

298 



238 



(C) 



ryrS S ryrJ AS 




2 hpf 
(64-cen) / ' 



2.25 hpf 
(12S-cell) 




Figure 3 Zebrafish ryr mRNA can be detected from very early developmental stages through to adulthood, (a) Temporal expression 
pattern of zebrafish ryrs using end-point PCR analysis and isoform specific primers. Total RNA was isolated from different developmental stages, 
as indicated, and reverse transcribed. 1.5 |jg of cDNA was synthesised and isoform specific primers for ryrla, ryrlb, ryr2a, ryr2b and ryrS used. 
Numbers in the right panel indicate the molecular mass of PCR products, p-actin was used as an internal control. The figure shows a 
representative result of replicates from at least three experiments, (b) Expression of zebrafish ryr2a and ryr2b in the brain (BR) and body (BO) of a 
3 months old adult zebrafish. p-actin was used as an internal control. Numbers in the right panel indicate the molecular mass of RT-PCR 
products, (c) Maternal expression of ryrS mRNA was observed as staining (arrows) within the dividing cells of 2-2.25 hpf embryos (64 and 128- 
cell stage). The specificity of the signal is compared to the ryrS sense probes. 



After the lateral migration of the slow muscle precursors 
is completed, the expression of ryrla mRNA is clearly 
visible in the superficial layer of slow muscle (Figure 
4b). There is conflicting evidence that by 24 hpf ryrla 
mRNA is confined either to slow muscle fibres [13], or 
can be detected at low levels throughout the somite [14]. 
In this study we found that ryrla mRNA co-localises 
with a known marker of slow muscle fibres at 24 hpf 



(Figure 4b). Furthermore, ryrla mRNA staining was not 
detected in the somites of smoothened (smo-/-) 
embryos, a mutant line which lacks slow muscle and the 
muscle pioneers (Figure 4c). Thus, we conclude that 
ryrla mRNA is expressed exclusively in embryonic slow 
muscle. 

We observed low levels of ryrlb expression in adaxial 
cells adjacent to the notochord beginning at 12-13 hpf 
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Figure 4 Zebrafish ryrla, ryrib and ryr3 mRNA is localised to embryonic skeletal muscle, (a) Expression of zebrafish myoD, ryrla, ryrlb 
and ryrS was examined using wiioie mount in situ iiybridisation. Expression of ryrla mRNA is detectable in tine adaxial cells of 1 1 hpf embryos 
(3-somite, arrows), whereas ryrlb mRNA expression is present in cells adjacent to the notochord from 12 hpf (6-somite stage) embryos. 
Expression of ryrS mRNA was only evident at 24 hpf, with the strongest staining observed in the anterior somites. The embryos at 1 1 and 12 hpf 
are orientated so that the anterior is to the left, (b) Cross-sections showing myoD, ryrla, ryrlb and ryrS in-situ mRNA hybridisation (above) and 
double immunostained with the F59 antibody/HRP labeling (below) in 24 hpf embryos, (c) ryrla, ryrlb and ryrS mRNA expression in wildtype 
(smo^^^), heterozygote (smo^^~) and homozygous (smo"^~) mutants at 24 hpf. There is no ryrla mRNA expression in the smo~^~ mutant. 



compared to wildtype and heterozygous embryos at 24 hpf. 
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(Figure 4a). At 18 and 24 hpf, strong ryrlb expression 
was detected throughout the somites in a pattern analo- 
gous to that of myoD, suggestive of its presence in both 
fast and slow muscle. Currently there is conflicting 
reports suggesting that ryrlb expression is confined 
either to the fast muscle [13] or appears in both fast 
and slow muscle [35]. In the later study the observation 
that ryrlb mRNA expression is found in both muscle 
types was attributed to an artefact arising as a conse- 
quence of contaminated tissue [35]. In order to clarify 
our findings we performed double-fluorescent labeling 
experiments, to observe ryrlb expression in the presence 
of a slow muscle marker. In contrast to previous work 
which used a ryrlb in situ probe that targeted the rela- 
tively conserved pore forming region at the C-terminus 
we generated a probe that recognised a more variable 
region at the N-terminus of the receptor [13]. Our data 
revealed that ryrlb mRNA is localised to both fast and 
slow developing skeletal muscle (Figure 4b, Additional 
file 1: Figure SI). 

Maternal expression of ryrS in the dividing cells of 2 
to 2.25 hpf (64- and 128-cell) embryos was confirmed 
by wholemount in situ hybridisation (WISH) (Figure 
3b); however, the low level expression of ryrS from 5.3 
to 18 hpf detected in the PGR analysis was not con- 
firmed by WISH (data not shown). In this study ryrS 
mRNA expression was detected in the skeletal muscle 
from 24 hpf using WISH (Figure 4a,b), with stronger 
expression in the anterior compared to developing pos- 
terior somites (Figure 4a). There is conflicting evidence 
that ryrS mRNA in zebrafish is expressed exclusively in 
slow and fast skeletal muscle from 14 hpf [11] or in 
many tissues including the CNS during the somitogen- 
esis [14,36]. The discrepancies in the expression patterns 
of ryrS mRNA may be explained in part by the target 
sequence used for probe synthesis in previous studies. 
The ryrS clone (019-D04-2) used previously targeted the 
conserved pore-forming region at the C-terminus of the 
receptor which shares high similarity to both the ryrla 
and ryrlb sequence, this raises the possibility that cross 
hybridisation with ryrl isoforms occurred [11]. In the 
current work the ryrS probe was designed to a more 
divergent region located in the N-terminus of the recep- 
tor. In this study ryrS mRNA expression was not 
detected in the Kuppfer's vesicles, adaxial cells at 12 hpf 
or in any region of the CNS at 24 hpf (data not shown). 
Double staining experiments revealed that ryrS mRNA 
expression was confined to the fast skeletal muscle at 24 
hpf (Figure 4b, Additional file 1: Figures SI and Addi- 
tional file 2: Figure S2). Thus we conclude that ryrS 
mRNA is expressed exclusively in the fast skeletal mus- 
cle of embryos at 24 hpf. 



ryr2a and ryr2b are expressed in the developing zebrafish 
CNS and heart, respectively 

Expression of ryr2a mRNA was detected initially in the 
CNS of zebrafish embryos starting at 24 hpf (Figure 5a). 
At 24 and 48 hpf, ryr2a mRNA expression can be 
observed at several sites in the developing forebrain (tel- 
encephalon and diencephalon), midbrain (mesencepha- 
lon) and hindbrain (rhombencephalon) of the zebrafish 
embryo (Figure 5a). A lateral view shows ryr2a mRNA 
expression in all seven hindbrain rhombomeres, rl-r7 by 
48 hpf (Figure 5b). The dorsal view clearly revealed two 
bilateral patches of ryr2a expression lying dorsal to the 
eyes which may correspond to the developing tegmental 
region (Figure 5c). Low level staining was also observed 
in the anterior portion of the spinal cord at 48 hpf (Fig- 
ure 5a). Analysis of ryr2a expression by PCR showed 
that 3 month old zebrafish display a strong ryr2b 
mRNA expression in the brain compared to the body 
(Figure 3b). However, the significant maternal expres- 
sion of ryr2a mRNA detected by PCR (Figure 3a) was 
not detected by WISH (data not shown). We conclude 
that zygotic ryr2a expression is expressed exclusively in 
the developing CNS of zebrafish embryos. 

The expression of ryr2b was compared to two well- 
characterised cardiac markers, nkx2.S and cmlc2 
[18,20,37] (Figure 6). Expression of nkx2.S, a homeodo- 
main transcription factor, is the earliest known marker 
of zebrafish precardiac mesoderm [20,37]. Cardiac myo- 
sin light chain 2 (cm/c2), another marker of myocardial 
precursor cells, is expressed throughout heart fusion 
and early heart tube assembly [18]. We show that 
nkx2.5 but not cmlc2 and ryr2b expression can be 
observed in the heart progenitor cells at 11 hpf; how- 
ever, by 14 hpf the bilateral expression of cmlc2 and 
ryr2b was apparent in the precardiac mesoderm. As 
reported previously [18], nkx2.S expression extends 
beyond the anterior tip of the notochord and the most 
posterior nkx2.5 cells {nkx2.5^) do not express cmlc2 
(cmlc2'). ryr2b expression appears to be expressed in the 
region corresponding to cmlc2^ consistent with its 
appearance in cells that will contribute to the myocar- 
dium. By 18 hpf the bilateral cardiac primordial cells 
expressing cmlc2 and ryr2b bend towards each other, 
make contact and begin to fuse (Figure 7a). The poster- 
ior portions fuse initially followed by anterior portions 
to create a central lumen and cardiac cone by 20 hpf. 
At 24 hpf, the heart has formed a linear tube with atrial 
and ventricular precursor cells and ryr2b expression 
appears to mirror that of cm/c2"^cells, suggestive of its 
expression throughout the heart at this stage. Using 
WICC, we observed RyR protein in the atrium and ven- 
tricle of the zebrafish heart at 48 hpf (Figure 7b). 
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Figure 5 Zebrafish ryr2a is expressed in the central nervous system. The spatial distribution of ryr2a mRNA in the zebrafish embryo was 
examined using whole mount in situ hybridisation, (a) Specific staining was observed initially in localised regions of the brain, specifically the 
telencephalon (T), diencephalon (D), mesencephalon (M), tegmentum (Tg) and rhomobomeres (rh) of whole mount embryos at 24 and 48 hpf 
(b) Lateral and (c) dorsal views (eyes removed) of the brain at 48 hpf revealed that ryr2a expression was present in each of the seven 
rhombomere segments, the tegmentum, trigeminal ganglia and anterior portion of the spinal cord. Embryos are orientated so that anterior is to 
the left. 



Discussion 

Structure of the zebrafish ryr genes and their products: 
comparison to other species 

Previous work identified a total of 14 genomic contigs 
for ryr from the zebrafish genome assembly [13]. Classi- 
fication of RyR sequences by radiation hybrid mapping 
suggested that there are at least five different zebrafish 
genes: ryr la, ryr lb, ryr2a, ryr2b and ryrS [13]. Our 
study has confirmed the sequences of the five zebrafish 



ryr genes, although currently the ryr2a gene annotation 
is still incomplete. Zebrafish, like other teleosts, have 
undergone a gene duplication event and appear to have 
retained two distinct copies of ryrl and ryr2. The teleost 
ryrl genes are differentially expressed in the skeletal 
muscle tissue and the receptors have distinctive Ca^^ 
binding sensitivities, suggesting that the genes have 
evolved to perform different physiological functions 
[13,38]. Here we report that the zebrafish genome 
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Figure 6 Zebrafish ryr2b is expressed in the differentiating cardiac muscle. Expression of zebrafish nkx2.5, cmlcl, and ryrlb was examined 
using wliole mount in situ liybridisation. Expression of nkx2.5 was observed in myocardial precursors from 1 1 lipf (3-somite stage) onwards, 
wliereas cmlcl and ryr2b were expressed in differentiating cardiac tissue from 14 lipf (10-somite stage) onwards. 11,12 and 14 lipf embryos are 
orientated so tliat tine anterior is to tine top, wliereas in 18, 24 and 48 lipf embryos anterior is to the left. 



contains two copies of the ryr2 gene, currently the only 
teleost in which this gene duplication can be observed. 
Our study has revealed that expression of ryr2a and 
ryr2h is confined to the embryonic nervous and cardiac 



tissue respectively. The differential expression of the 
ryr2 supports the idea that the gene products have dif- 
ferent roles within these tissues. The reported sub-func- 
tionalisation of the ryr2 in the zebrafish embryo in this 
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Figure 7 Zebrafish ryr2b mRNA is expressed during early heart formation and the RyR protein throughout the heart at later stages. 

(a) Expression of nkx2.5, cmlc2, and ryr2b in 18 lipf (18-somite stage) was observed in tine bilateral cardiac primordial cells. At this stage the 
bilateral cardiac primordial cells make contact and begin to fuse (arrows). The posterior portions fuse initially followed by anterior portions to 
create a central lumen and cardiac cone by 20 hpf (22-somite stage). A linear heart tube has formed by 24 hpf and ryr2b and cmlc2^ cells are 
expressed throughout the heart at this stage. 18- and 22-somite stage embryos, in which the tail was removed but the yolk sac left intact, were 
orientated so that anterior is to the top. Dorsal view of flat mounted embryos at 24 hpf. (bSfe) Brightfield images of the anterior portion of the 
zebrafish embryo with the yolk sac still intact and dorsal side uppermost. Embryos were fixed and stained at 48 hpf to reveal either myosin (c-d) 
or RyRs (f-g), the images within the white boxes are shown in greater detail (d & g). By 48 hpf RyRs are expressed throughout the two 
chambers, atrium (A) and ventricle (V) of the heart. Images were taken using XI 0 (b,c & e,f) and X20 (d & g) objectives. 
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study will facilitate the study of this receptor in the 
development of nervous and cardiac tissue. 

Our work has revealed that the primary sequence of 
the zebrafish RyRs contains many of the conserved 
domains associated with the regulation and function of 
this intracellular ion channel in other species, mostly 
notably humans. The primary sequence of the Ca^^-con- 
ducting pore domain was found to be extensively con- 
served between zebrafish and other species examined. 
Our study of RyR primary structure supports previous 
work to show that the biophysical properties and phar- 
macological regulation of the zebrafish RyRl is similar 
to its mammalian homolog [39]. Taken together this 
data reveals that the zebrafish RyR functions in a similar 
manner to those found in mammals, this information is 
of significant relevance to work using the zebrafish as a 
model for human disease. However, differences in the 
predicted primary sequence of the zebrafish ryanodine 
receptor proteins and those of other species were also 
recorded. There are reports species-specific differences 
in the Ca^^ regulation and single channel conductance 
of the RyRl channel [39]. Further comparative analysis 
of the RyR family will provide a better insight into the 
physiological functions of the receptor at a tissue or 
whole organism level. 

The expression of ryr genes during early development, 
up to and including axis formation 

Calcium signalling is required throughout development; 
however the signal pathways have not been well defined. 
The prevailing view is that IPsR-driven Ca^^ signals 
have a major role in axis formation prior to organogen- 
esis whereas RyR-induced Ca^"^ signals are necessary for 
the later aspects of tissue specific differentiation (e.g. 
muscle formation) [3-5,7,9,40]. Our current study 
revealed strong maternal expression of ryr2b and ryrS 
genes during cleavage, with low levels of ryrS expression 
detectable throughout the blastula and gastrula periods. 
The significance of early ryr mRNA expression remains 
to be determined; however, it raises the possibility that 
RyR-generated Ca^^ signals act in development prior to 
10 hpf. The characterisation of expression in zebrafish 
establishes a basis for future experimental work aimed 
at determining the action of RyR-induced Ca^"^ signal- 
ling events in early embryonic patterning. 

The expression of ryr genes in developing skeletal muscle 

Several studies have implicated RyR function in organo- 
genesis, particularly in striated muscle development (as 
reviewed in [9]). Our study has used both mutant lines 
and double staining to establish that RyRs are differen- 
tially expressed in slow (ryrla), fast (ryrS) and both types 
{ryr lb) of developing skeletal muscle during the segmen- 
tation period. In E-C coupling within mammalian skeletal 



muscle, RyRl is directly coupled to a voltage gated-Ca ^ 
channel (VOC) on the sarcolemma. Activation of VOCs 
via membrane depolarisation then triggers the opening of 
the ryanodine receptor (RyRl) and release of Ca^"^ from 
the sarcoplasmic reticulum (SR) stores. Zebrafish have 
non-Ca^^-conducting voltage-gated Ca channels that 
have evolved solely as voltage sensors to trigger opening 
of the RyR [35]. The central subunit of the VOC 
(Cavi lais) ^cts as the pore, selectivity filter and voltage 
sensor. In zebrafish two Cayi i^is genes {zf-als-3. and zf- 
als-b) have been identified. The Cayi.iais gene products 
are proposed to interact in a tissue specific manner with 
the ryrl genes; that is zf-als-a and ryr la are expressed in 
slow muscle whilst zf-als-b and ryr lb are confined to 
fast muscle. However, our data suggest that the situation 
is not quite as clear cut as first proposed because ryr lb is 
not expressed exclusively in fast muscle but is also 
located in slow muscle. In addition, Ca^^ release can also 
be regulated by the ryrS gene product which is also 
located in the fast skeletal muscle. RyR3 is proposed to 
act as an uncoupled calcium-induced calcium release 
(OCR) channel to propagate the Cd?^ signal [41]. There- 
fore we propose that in the developing fast muscle zf- 
als-b and ryr lb act together to generate an increase in 
[Ca^"^]i, with the ryrS gene product acting to amplify the 
signal. Our data in the embryonic fast muscle showed 
that ryr lb expression occurs prior to ryrS and suggests 
that the RyR-generated [Ca^"^]i increase occurs initially 
via ryr lb with the proposed amplification step via ryrS 
developing subsequently. In zebrafish the role of the RyR 
in E-C coupling within the developing slow muscle 
appears more complex. Our data reveals that both ryrla 
and ryr lb are expressed in the developing slow muscle, 
but ryrS is not. This presents the possibility that zf-als-a 
could couple to both ryrla and ryr lb and raises the issue 
of whether amplification of the Ca^^ signals occurs in 
this tissue and, if so, how is this achieved. Clearly there is 
still much to understand about the maturation of depo- 
larization-induced Ca^"^ signaling and its role during ske- 
letal muscle differentiation in vivo. 

The expression of ryr genes in the developing nervous 
and cardiovascular systems 

Our data has revealed that there are two ryr2 genes, 
ryr2a and ryr2b, which are exclusively expressed in 
either the developing nervous system or cardiac tissue, 
respectively. Studies in mammals revealed that RyR are 
expressed in the developing brain and that RyR- 
mediated Ca^^ signals may have a role in neuronal dif- 
ferentiation and neurite outgrowth [42-44]. All three ryr 
genes are expressed within the embryonic mouse brain; 
however, from postnatal day 7 onwards ryr2 becomes in 
the major isoform [45]. The postnatal changes in RyR 
expression in mouse brain correlate with a period of 
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neuronal differentiation and may therefore be important 
in establishing [Ca^'^]i homeostasis in maturing neurons. 
In zebrafish ryr2a expression is localized to specific 
regions of the developing brains. In these regions ryr2a 
is likely to regulate neuronal Ca^"^ signaling and there- 
fore play a role in CNS development. 

In mature cardiac muscle, Ca^"^ signals are generated 
by CICR via the activation of VOCs and the cardiac RyR 
(RyR2). Knockout mice which do not express ryr2 initi- 
ally display spontaneous rhythmic contractions of the 
heart at embryonic day 9 (E9) but no heart beat by day 
10 (ElO) [4]. Furthermore, RyR-mediated Ca^"^ release 
does not play a significant role in the [Ca^"^]i changes 
observed within the heart of new born rats [46]. Thus in 
mammals it appears that the RyR2 does not contribute 
to the onset of contractile activity at very early embryo- 
nic stages, but is important for the subsequent matura- 
tion and development of the heart in vivo. The zebrafish 
cardiac ryr gene {ryr2b) is expressed exclusively in the 
developing heart tissue (precardiac mesoderm) from 14 
hpf, 8 hours prior to the onset of cardiac contraction at 
22 hpf, and may well contribute to early cardiac devel- 
opment. Investigation of ryr2 function during mamma- 
lian development is complicated by the fact that a single 
ryr2 gene is expressed in several tissues [1]. The sub- 
functionalisation of the ryr2a and ryr2b genes in the 
zebrafish embryo provides an excellent system to study 
individual receptor function in neuronal and cardiac tis- 
sues during vertebrate development. 

Conclusions 

This study has provided a comprehensive overview of 
the spatial and temporal expression of the ryr gene 
family in developing zebrafish embryos. This family of 
Ca^"^-release channels are expressed predominantly in 
developing skeletal, cardiac and neuronal tissue, suppor- 
tive of the view that RyRs function is relevant to later 
development events, such as tissue differentiation. In 
addition, the study has also revealed that maternal ryr 
mRNA is present in the very early embryo, suggestive of 
a function for this receptor prior to organogenesis. Rya- 
nodine receptors have been implicated in human disease 
and the zebrafish is an important vertebrate develop- 
mental model which will facilitate work in this area. 
Future work will explore the function of RyR-regulated 
Ca^"^ signal pathways during zebrafish embryogenesis. 

Additional material 



Additional file 1: Figure S1. ryr1a, ryr1b and ryr3 mRNA is 
differentially expressed in the developing skeletal muscle. Double 
labelling in wholemount embryos at 24 hpf was performed using probes 
to myhzl, ryrla, ryrlb and ryr3 and Fast Red as a substrate (red) followed 
by immunostaining using the F59 antibody and a fluorescent secondary 
(green). Images show Z-stacks of whole-mount double-labelled embryos 



(top row) or sections (bottom row) from dissected embryos. Cross- 
sectional images revealed that ryrla co-localised exclusively with slow 
muscle staining, whereas myhzl, a marker of fast muscle, was not 
expressed in the slow muscle. Furthermore ryrlb expression could be 
observed throughout both muscle types whereas ryr3 did not co-localise 
with the slow muscle staining and appeared to be expressed exclusively 
in the fast muscle. Scale bars = 20 pm, unless otherwise indicated. 

Additional file 2: Figure S2. ryr3 mRNA expression is confined to the 
fast muscle fibres throughout the myotome at 24 hpf. Double 
fluorescent cross-sections (1-3) were prepared by labelling a 24 hpf 
zebrafish embryo with a fluorescence substrate (i.e. FastRed; red) for ryrS 
in situ hybridisation and the F59 antibody (green) for immunostaining. 
The position of cross-sections 1-3 are illustrated in the 24 hpf ryr3 WISH 
labelled embryo (top left) which has been stained with BM purple and 
laterally orientated with anterior to the left. Scale bars = 20 ^im, unless 
otherwise indicated. 
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